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Abstract 

We calculate the inclusive decays B — > Xgif and B — > XsPc using factorization 
assumption. To investigate the bound state effect of the decaying B meson in 
these inclusive decays we take into account the motion of the b quark using 
a Ganssian momentum distribution model. The resulting correction to free 
quark decay approximation is aronnd %6 at most. Utilizing a potential model 
evalnation of the ratio of the decay constants /^^//^, it is shown that the ratio 
R = T{B —> As?7c)/r(B ^ Xgf^) can be used as a possible test of factorization 
assumption. 
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Exclusive and inclusive nonleptonic B decays to charmonium states are of special interest 
theoretically and experimentally. These decay channels, among other things, provide a 
powerful testing ground for color suppression and factorization in hadronic B decays [|^. 
At the same time, exclusive modes of two body B decay into K meson resonances and 
charmonium states can provide an alternative examination of models for the treatment of 
hadronic form factors |@] . 

In this work, we focus on the inclusive two body decays B —>• Xs'ijj and B —> X^r^c, 
where Xg is a final state hadron containing a strange quark. There is no experimental data 
on the latter decay at this time. However, as we point out, the eventual measurement of 
this inclusive decay channel can be used to test the validity of the factorization assumption 
in nonleptonic B decays. In fact, one can show that the ratio R = T{B Xgrjc)/T{B 
Xg'ip), calculated by using factorization, is independent of the QCD corrections and the scale 
ambiguity of the Wilson coefficients. In this context, R depends only on the ratio of the 
decay constants fr^c/ftp for which we use an improved estimate obtained in a previous work 


The inclusive decays B —>• Xg'ip^ric) are usually approximated with the free quark decays 
b —> s'ijj{ric). To improve upon, in the present paper, we estimate the correction to this 
approximation by taking into account the motion of the b quark inside the B meson. For 
this purpose, we use a one parameter Gaussian momentum distribution for the b quark which 
has previously been applied to inclusive semileptonic @] , rare dileptonic 0 and nonleptonic 
B decays [^(commonly known as ACCMM model in the literature). We present results for 
a range of the model parameter obtained from fits to experimental data. 

Neglacting penguin operators, the relevant effective Hamiltonian for B Xs'ip{r]c) can 
be written as: 


fi.ff = 


Gi(/r)cV(l - - j5)c’ 

+G2(p)sV(1 - - 75)<^ 


+ H.C. , 


( 1 ) 


where i and j are color indices and C'i(p) and C' 2 (/r) are QCD improved Wilson coefficients. 
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One then can use a Fierz transformation to write (1) in the following form: 




[(C2(h) + iC'i(h)) sV(l - 75 )&*c^7m(1 - 75)c^' (2) 

+ Oi(/i)sV(l - + H.C. , 

where T“ (a = 1..8) are generators of SU{3)co\or- We note that the hrst term in (2) is 
the product of two color singlet currents but the second term consists of two color octet 
currents. In the factorization assumption, only the color singlet quark current contributes 
to the production of colorless cc hnal state -0 and rjc- Using the dehnition for the decay 
constant {f^) for the vector meson 'ip\ 


=< 0|c7^c|i/> > , 


(3) 


(e^ is the polarization vector of ip) the effective Hamiltonian for B —>■ Xgip is obtained as 
follows: 


= cusY{i - . (4) 

where 

C = '^vxv,. (cm + ^cm) ■ ( 6 ) 

Similarly, utilizing the dehnition for the decay constant for the psuedoscalar meson rjc: 

fvcQ^ =< 0|c7'^75c|7c(g) > , (6) 

results in the following effective Hamiltonian for B —>■ XsTjc decay: 

= (7) 

The Wilson coefficients Ci and C 2 in (5) are calculated to the next-to-leading order in 
reference |]^ resulting in 02 = C' 2 (/i) -|- l/3C'i(/i) = 0.155 for /x = m;, ~ 5 GeV. However, the 
branching ratio for the exclusive decay B —>■ Kip obtained from (4) requires 02 to be roughly 
by a factor of two larger than the above value in order to agree with the experimental data 


BR{B+ K+iP) = (0.101 ± 0.014)% . 
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This discrepancy between theoretical prediction and measurement could be due to two fac¬ 
tors. On one hand, the /i-dependence of the Wilson coefficients which arises from short 
distance QCD results in a signihcant uncertainty in the calculated decay rate in the con¬ 
text of factorization. In fact, phenomenologically, 02 is treated as a free parameter to be 
determined from experiment |^. On the other hand, one may question the validity of the 
factorization assumption which allows to infer eq. (4) from the effective Hamiltonian (2). 
In other words, the second term in (2) which is nonfactorizable could have a signihcant 
contribution to the matrix element |]^.. To disentangle these two factors and examine the 
factorization assumption, the ratio of the inclusive decays R can serve as a crucial test¬ 
ing ground. Aside from the cancellation of the Wilson coefficients in R, this ratio is also 
free from the nonperturbative hadronic uncertainties which is usually associated with the 
theoretical calculations of exclusive decays . 


Using (4) and (7), one can calculate the decay rates r(6 —>• s'lp^rjc))'. 

C^fl 


T{b s'iIj) = 






X 


2/’2i 2^ 2/02 2^1 4 0‘ 

mb[mb + m^) — m^[2mb — m^) +mg — 2m. 


( 8 ) 


T{b^ sr]c) = 


where 


\{ml - m^f - + m^) 


Sirmb' 


9 {x,y,z) = 


2 2 2 






2 2 n 

, y 

x"^ 


(9) 


( 10 ) 


and mb and mg are b and s quark masses, respectively. The inclusive decay rates r(i? — 
^s'll’iVc)) are usually approximated by eqs. (8) and (9). However, in this work we estimate 
the bound state corrections to this approximation by taking into accout the motion of the 
heavy b quark inside the B meson. We follow the ACCMM method ^ which incorporates 
the bound state effect in semileptonic B decays by assuming a virtual b quark inside B meson 
accomponied by an on-shell light quark. In the meson rest frame, the energy-momentum 
conservation leads to the following relation for b quark mass W: 
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(11) 


iy2(p) = ^2 _ 2mB^p‘^ + m2 , 


where is the light quark mass and p is the 3-momentum of the b quark. Following 
reference IQ, we also consider a Gaussian momentum distribution for the Fermi motion of b 
quark: 

^(P) = . (12) 

The model parameter pp determines the distribution width, and is related to the average 
momentum < p >. 

At this point we would like to remark on the consistency of the above model with heavy 
quark expansion. Let us consider the average b quark mass dehned as: 


rPmax 

rni= W(p)if>(p)p dp , 


( 13 ) 


where Pmax is the maximum kinematically allowed momentum. Using eq. (13), one can 
derive an expansion of the B meson mass in powers of as follows: 


_ , 3p^ / 1 

rriB = mb-\ -^ U 


TT Arrih 


m 


2 ; ’ 


(14) 


in which rUg = 0 is assumed. A comparison of eq. (14) with the usual heavy quark expansion 
formula for heavy-light mesons, i.e. 


T Xl+dM^2 I 1 \ 

= + + , ( 15 ) 

reveals that once rub is identihed with the mass of the heavy quark tuq in eq. (15), the 
nonperturbative parameters A and Ai of the heavy quark expansion and the model parameter 
Pf are connected as follows: 


- ^ JA 

sft 2 


( 16 ) 


The ACCMM model does not provide a corresponding term for the nonperturbative param¬ 
eter A 2 (c^M = 3,-1 for psuedoscalar and vector mesons, respectively) which is due to the 
spin interactions. This could be considered as a shortcoming of the model. However, the 
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constraint imposed by eq. (16), i.e. Ai = —Svr/SA^, is in reasonable agreement with qnoted 
valnes for these parameters 0- It is in this sense that we consider the above model to be 
consistent with heavy qnark symmetries. 


To incorporate the effects of the motion of the b qnark in the inclnsive decays B —>■ 
Xs'ipirjc), we replace the h qnark mass mb in eqs. (8) and (9) with W{j>) dehned in eq. (11) 
and integrate over the kinematically allowed range of the b qnark momentnm p, i.e. 


fPmax 

T{B Xs^p{r|c)) = T{b^ s^p{r]c))m,=w(p)(|>{p)p dp . 
0 


( 17 ) 


As a resnlt, the sensitivity to the heavy qnark mass mb is replaced with the model parameter 
Pf and the light qnark mass dependence. There are varions determinations oi pp from 
hts to semileptonic B decays and also from comparison with the heavy qnark effective 
theory approach. In reference 0, A = 0.55 ± 0.05 GeV and Ai = —0.35 ± 0.05 GeV^ have 
been extracted from GLEO data on inclnsive semileptonic B Xiui decay. A comparison 


with eq. (16) leads to pp ~ 0.5 GeV. However, in reference [^], a smaller central valne 
Pf = 0.271q27 is obtained by an AGGMM model analysis of the ARGUS resnlts for the 
lepton energy spectrnm of R —> X^iPe decay. In order to investigate the sensitivity of onr 
estimates to the model parameter pp, we present resnlts for pp = 0.3 and pp = 0.5 GeV. 
Using eq. (14), these valnes of pp correspond to nib = 4.92 GeV and nib = 4.65 GeV, 
respectively. 

Following the prescription of eq. (17) and introdncing the notation 


T {B ^ Xsi^iPc)) = 


f2 

8tt 


N, 


PiVc) ) 


( 18 ) 


we obtain = 7.64(8.14) GeV and = 44.29(48.18) GeV^ for pp = 0.3 GeV (the hrst 
nnmber is calcnlated by taking the strange qnark mass = 0.55 GeV and the second 
nnmber in parenthesis is resnlted from nsing = 0.15 GeV). Gomparing these resnlts with 
the case where mb = nib = 4.92 GeV is inserted in eqs. (8) and (9), i.e. = 7.55(8.04) GeV 
and Nri^ = 43.55(47.42) GeV^, indicates that only a small correction of order %1 — 2 arises 
from considering this bonnd state effect. On the other hand, a larger valne of the model 
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parameter pp = 0.5 GeV (which is compatible with heavy quark effective theory approach) 
results in = 5.85(6.35) GeV and = 33.14(36.77) GeV^. A comparison with the decay 
rates obtained from eqs. (8) and (9) by using mb = fUb = 4.65 GeV, i.e. = 5.57(6.06) 
GeV and Njj^ = 31.01(34.62) GeV^, reveals a larger bound state corrections of order %5 —6. 
We also notice that the corrections are almost independent of s quark mass m^. 

As we observe from eq. (18), the ratio R = r{B —> XgPc)/^{B —> Xgip) is independent 
of C in the context of factorization assumption, i.e. 


p _ XsTjc) _ ^Vcfvc _ /c o n 1 


(19) 


where our error estimate in the numerical factor represents the variation of the model pa¬ 
rameter pp in the range 0.3 to 0.5 GeV and also the s quark mass mg from 0.15 to 0.55 
GeV. 

To evaluate the ratio of the decay constants in (19), we use the potential model relations 
which relate these form factors to the value of the meson wavefunction at the origin: 


U = 


12 


m. 


^17.(0) , 


r?c 


( 20 ) 


The common assumption in the literature is that the wavefunction of the psuedoscalar and 
vector mesons at the origin are more or less identical. However, in reference 0, based on a 
simple perturbation theory argument, this ratio was estimated to be 

1^. 




= 1.4 ±0.1 . 


Gonsequently, the ratio of the decay constants in (19) can be written as: 


( 21 ) 


J Vc 


1^ 


Vc 


fp rnr,^m^ |^p(0)|^ 

Inserting (22) into (19), we obtain: 

TjB^ Xgp,) 
T{B ^ Xg-ilj) 


0.15 ±0.01 GeV 


-2 


( 22 ) 


= 0.87 ±0.06 


(23) 
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We would like to point out again that the ratio R is free of hadronic model uncertain¬ 
ties which are normally encountered in calculating exclusive decays. The inclusive decay 
branching ratio BR{B —>■ Xs'ip) = (1 ± 0.25) x 10“^ has been extracted from experimental 
data for direct production of "0 in nonleptonic B decays [|^. Therefore, a measurement 
of the inclusive B —> Xgrjc decay along with reducing the error bar in the experimental 
result for B —>■ Xg'ip can serve as an alternative test of the factorization assumption which is 
used in deriving eq. (23). A signihcant deviation of the experimental value for R from the 
theoretical prediction in eq. (23) could be an indication of the failure of the factorization 
assumption. We would like to emphasize that the cancellation of the Wilson coefficients in 
the ratio R results in a signihcant reduction of the uncertainty due to QCD corrections and 
scale dependence. 

In conclusion, using factorization, we calculated the ratio of the inclusive nonleptonic B 
decays to a hadron containing a strange quark plus psuedoscalar and vector cc mesons. The 
bound state effect due to the motion of the b quark inside the B meson was also considered. 
We used an improved estimate of the ratio of the dacay constants for rjc and 'ip was obtained 
based on potential model and pertubation theory argument. Once the experimental results 
OYi B ^ XgTjc decay are available, a comparison to the theoretical prediction presented 
in this paper could serve as a test of factorization assumption. Finally, we would like to 
emphasize that the pattern of deviations from factorization can give us important clues on 
the nonfactorizable contributions, and therefore, it is important to examine the departures 
from this approximation using various experimental data. 
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